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Abstract. We present a 3d code for the dynamical evolution of a multi-phase 
interstellar medium (ISM) coupled to stars via star formation (SF) and feedback 
processes. The multi-phase ISM consists of clouds (sticky particles) and diffuse gas 
(SPH): exchange of matter, energy and momentum is achieved by drag (due to 
^— ' ram pressure) and condensation or evaporation processes. The cycle of matter is 

, completed by SF and feedback by SNe and PNe. A SF scheme based on a variable 

' SF efhciency as proposed by Elmegreen & Efremov (1997) is presented. For a Milky 

Way type galaxy we get a SF rate of ~1M0 yr~^ with an average SF efficiency of 
-5%. 
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^ ! 1. Introduction 

o 

CN I So far in 3d models of galaxies the ISM is mostly described either 

as a diffuse phase by SPH (e.g. Hernquist & Katz, 1989) or as a 
^ [ clumpy phase by sticky particles (e.g. Theis & Hensler, 1993). Alter- 

natively, chemo-dynamical models use a multi-phase ISM, but were 
Q^. usually restricted due to spherical or axisymmetric systems (e.g. Theis 

^ I et al., 1992; Samland et al., 1997). In order to extend the chemo- 

^ ■ dynamical models to three dimensions we combine both treatments 

"to . in a particle based code: the hot diffuse gas phase is described by a 

I SPH formalism, whereas the cold molecular clouds are represented by 

sticky particles. The coupling between the gaseous phases is achieved 
^ I by condensation/evaporation (C/E) and by a drag force due to ram 

^ • pressure. Energy is dissipated by cloud-cloud-collisions or by radia- 

tive cooling. Furthermore, stars are formed in clouds and the stars 
return mass and energy to the gas by feedback processes (supernovae, 
planetary nebulae). 

In SPH codes the SF is generally based on the Schmidt law, i.e. 
the SF rate depends on gas density and a characteristic time scale. 
For sticky particles the SF can be coupled to cloud-cloud collisions, 
or a single cloud forms stars with a constant SF efficiency. Here, a 
SF scheme using a different approach is presented: the process of SF 
is treated individually for each cloud. The SF efficiency is dependent 
on local properties of the ISM and the star forming clouds, thereby 
enabling self-regulation of SF by feedback. 

© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 



harfst.tex; 2/02/2008; 2:47; p.l 



2 



S. Harfst et al. 




Figure 1. The SF scheme (see text for details) 



2. Numerical Treatment 

The gravitational force for all particles is calculated with a TREE 
algorithm. A dark matter halo can be modeled as an external, static 
potential (static halo) or by means of particles (life halo). The inte- 
gration of the system is done with a leap-frog scheme. Individual time 
steps are used for each particle. 

The drag force and the mass exchange rates for C/E are calculated 
following Cowie et al. (1981) using a mass-radius-relation for the clouds 
based on observations (e.g. Rivolo & Solomon, 1988). Effects of both 
processes are individually determined for each cloud using the local 
density, temperature and velocity of the hot gas. Energy can be dissi- 
pated by inelastic cloud-cloud collisions (Theis & Hensler, 1993) and 
by radiative cooling (Bohringcr & Hensler, 1989). 

Clouds can transform some of their mass into new stars. In Fig. 1 
this process is drawn schematically. Each cloud is inactive for a given 
time Tia (to). An embedded star cluster is then formed (ti) with the SF 
efficiency e being a function of cloud mass and gas pressure according 
to Elmegreen & Efremov (1997). The energy released by SNe is calcu- 
lated using a multi-power law IMF (Kroupa et al., 1993). This energy 
injection disrupts the cloud into smaller fragments (t2). The time for 
disruption T^is {^2 — ti in Fig. 1) is determined from the energy input: 
a self-similar solution is used to calculate the SN shell expansion. It is 
assumed that the cloud disrupts when the shell radius equals the cloud 
radius. The fragments are then placed symmetrically on the shell with 
velocities corresponding to the expansion velocity at T^is- The energy 
not used for cloud disruption and the mass ejected by SNe is returned 
to the hot gas phase (SPH particles). Additionally, the mass returned 
by PNe is added to surrounding cloud particles. 
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Figure 2. In the left plot the temporal evolution (solid line) and the average (dashed 
line) of the SF rate are shown for a Milky Way type galaxy. The right plot shows the 
average SF efficiency (solid line) and the range of individual SF efficiencies (shaded 
area) with time for the same model. 



3. The model and results 



Initially a galaxy similar to the Milky Way was set up using the galaxy 
building package described by Kuijken & Dubinski (1995). A cloudy 
gas phase is added by treating every tenth particle from the stellar 
disk cloud. The total mass in clouds is Mci,tot « 3.4 • 10^ M© and 
each cloud is assigned a time of inactivity Tja between and 200 Myr. 
Finally, a hot gas halo with a total mass of Mhot,tot ~ 2 • 10^ Mq is 
added, which is designed to be in hydrostatic equilibrium in the halo 
potential. 

In Fig. 2a the resulting SF rate is shown, which is on average 
^IMQyr"-*^. The SF rate is slowly decreasing with time due to 1) the 
decreasing mass in the cloudy phase of the ISM and 2) the decreasing 
average SF efficiency (Fig. 2b). The decrease of e is correlated with a 
decreasing M^ as expected from Elmegreen &; Efremov (1997). After 
700 Myr the average SF efficiency is less than 5% and the maximum 
efficiency is always lower than 20%. From such a low SF efficiency we 
do not expect any non-induced globular cluster formation in the disk. 

In Fig. 3a the mass return rate due to stellar death is shown: the 
SNe rate is about 10% of the SF rate and in agreement with the mass 
fraction of massive stars. The mass return rate by PNc is approximately 
equal to the SF rate due to a contribution of older stars in the initial 
model. 

Condensation and evaporation rates are shown in Fig 3b. The evapo- 
ration rate is about 10 times higher than the condensation rate because 
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Figure 3. Mass flow rates. The left plot shows the temporal evolution of the mass 
return rate by PNc (solid) and SNc (dashed). In the right plot condensation (dashed) 
and evaporation (solid) rates are shown. 



the temperature of the hot gas amounts to several 10^ K. While the SF 
establishes an equilibrium on a rather short time scale C/E have not 
reached an equilibrium state after 700 Myr. This is in agreement with 
the analytic work of Koppen et al. (1998). They predict an equilibrium 
state only after a few condensation time scales of evolution, which is a 
few Gyr in this model. 
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